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Abstract 

We calculate the elliptic flow of charged particles in Pb-Pb collisions at ^/s = 2.76 
TeV in relativistic viscous hydrodynamics. The recent data of the ALICE Collab- 
oration on the elliptic flow as function of the centrality can be very well described 
using the hydrodynamic expansion of a fluid with a small shear viscosity rj/s = 0.08. 
The elliptic flow as function of the transverse momentum shows systematic devia- 
tions from a hydrodynamic behavior in the small momenta region p± < SOOMeV. 
It indicates that a non-negligible contribution of non-thermalized particles from jet 
fragmentation is present. 

PACS: 25.75.-q, 25.75.Dw, 21.65.Qr 

The elliptic flow measurements in heavy-ion collisions at the Large Hadron 
Collider (LHC) have been presented fTJ. Generally, the results are similar as 
observed at lower energies at the Relativistic Heavy Ion Collider (RHIC). 
The elliptic flow coefficient V2{p±) as function of the transverse momentum 
is increasing with p± and saturates at higher momenta. The elliptic flow as 
function of the centrality of the collision reflects the initial eccentricity of 
the fireball at each centrality. The elliptic fiow is generated in a collective 
expansion of the dense matter created in the collision, the comparison to 
model calculations can provide valuable information on the dynamics of the 
collisions, on the equation of state and the transport coefficients of the dense, 
hot matter [2J3]. 

The integrated elliptic flow as function of the collision centrality can be used 
to extract the properties of the fluid in the expanding flreball. Analyzes of 
the elliptic flow in heavy-ion collisions at RHIC energies indicate that the 
dense matter is an almost perfect fluid [1|5] . The ALICE data from Pb-Pb 
collisions at ^/s = 2.76TeV give smaller V2 than expected from perfect fluid 
hydrodynamics with hadronic rescattering [6j. The comparison of viscous hy- 
drodynamic results with data at RHIC and LHC show no noticeable change 
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in the shear viscosity coefficient with the energy [7]. Analysis of the LHC 
data at different centrahties in terms of the Knudsen number points towards a 
small value of the shear viscosity [8j. The observed elliptic flow can be partly 
reproduced in a hadron rescattering model with a short formation time [S]. 
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Fig. 1. Charged particle multiplicity density in pseudorapidity per participant pair 
for Pb-Pb collisions at y/s = 2.76 TeV from the Glauber model (dashed line), and 
from viscous hydrodynamic calculations (solid line) together with the ALICE data 

m- 

In the following we present an analysis of the ALICE results for the charged 
particle elliptic flow in the second order viscous hydrodynamic model with 
shear and bulk viscosities [TT] . Particle production in the central region at the 
LHC is described in the 2+1-dimensional, boost-invariant hydrodynamics. The 
shear viscosity to entropy ratio is rjjs = l/in and the bulk viscosity is non- 
zero in the hadronic phase with (/s = 0.04. The initial entropy density of the 
fireball in the transverse plane is taken from the Glauber model, at the impact 
parameter b 



s{x,y) = So 



[1 - a)ppart{x, y, b) + 2apbi„{x, y, h) 
[l-«)Ppart(0,0,0) + 2apMn(0,0,0) 



with ppart and pun the participant nucleons and binary collisions densities. The 
details of the hydrodynamic model calculation and the parameters for Pb-Pb 
collisions at the LHC can be found in p^fTT] . The statistical emission at the 
freeze-out and resonance decays are done using the THERMINATOR event 
generator [13]. The parameter a = 0.15 reproduces the centrality dependence 
of the density of charged particles in pseudorapidity ^ measured by the 



2 



0.4 p 



0.35 r 



: Pb-Pb 2.76TeV ALICE Data 



0.25 r 



0.3 '- 



10-20% 
20-30% 
40-50% 



0.2 '- 




0.15 r 



0.1 '- 



0.05 h 




0.5 



1 



1.5 



2 



2.5 3 
[GeV] 



3 



Fig. 2. Elliptic flow coefficient of charged particles as function of transverse momen- 
tum measured by the ALICE Collaboration ^ (symbols) together with the results 
of viscous hydrodynamics (lines) for three centrality classes. 
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(dashed line in Fig[T]), where dNpp/dr] is the multiplicity density in inelastic 
proton-proton collisions, Npart and Nf,in are the average numbers of partici- 
pants and of binary collisions for a given centrality. The entropy density at 
the center of the fireball sq in Eq. [T]is adjusted to reproduce the multiplicity 
in the most central Pb-Pb collisions, taking into account the entropy produc- 
tion in the viscous hydrodynamic evolution. The results of the hydrodynamic 
calculation (solid line in Fig. [1]) follow closely the Glauber model input and 
reproduce the data on the centrality dependence of the multiplicity [TOj . 

For each centrality class the elliptic flow of charged particles as function of 
the transverse momentum is calculated in the reaction plane (Fig. [2]). For 
central collisions (10-20%) the results of the hydrodynamic calculation follow 
the data up to p± ~ 2GeV. For more peripheral collisions (40-50%) the data 
show a saturation of the increase of V2 with p± around 1.2GeV, while the 
calculated V2{p±) keeps growing with the transverse momentum. The change in 
the behavior of the elliptic flow coefficient with pj_ is expected. The production 
of particles with small momenta pj_ < l-2GeV is dominated by the statistical 
emission from a thermalized fluid. The asymmetry of the flow of the fluid 
comes from the hydrodynamic evolution of a source with an initial azimuthal 
eccentricity. On the other hand, particles at high momenta originate from jet 
fragmentation Differences in the path length in different directions in 
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the fireball, lead via jet quenching to an azimuthal asymmetry of the final 
particles. While it is expected that the hydrodynamic results deviate from the 
data at high transverse momenta, we observe also a systematic deviation of 
the data points from the calculation in the range 200MeV< p± < SOOMeV. 
We discuss about this effect latter. 

In Fig. |3] is shown the integrated elliptic flow f 2 for different centralities of 
the collision (|?7| < 0.8, 200MeV< p± < 5TeV). The hydrodynamic model 
using Glauber model initial sources, with shear viscosity i]/s = 0.08 and bulk 
viscosity (/s = 0.04 reproduces the experimental observations very well. It 
means that, for the chosen initial eccentricity, the expansion of the fluid with 
a minimal shear viscosity describes the data. The analogous calculation for 
Au-Au collisions at 200GeV describes well the RHIC data for the integrated 
V2 using the same value rj/s = 0.08. The precise value of the shear viscosity 
coefficient compatible with the data depends on the model of the initial eccen- 
tricity, with color glass condensate initial conditions leading to larger values 
oi rj/s imS]. We note that the Glauber binary collisions profiles of the initial 
energy density used in Ref. result in up to 20% larger initial eccentricities 
in central collisions than Eq. [H which would imply a larger value of rj/s to 
reproduce the data. 

Despite the apparent success of the hydrodynamic description of the average 
elliptic flow in Pb-Pb collisions at the LHC as presented in Fig. [3l the result 
must be taken with caution. The differential elliptic flow coefficient V2{p±) is 
not described by the hydrodynamic model, neither at high momenta, which 
is natural, nor at small momenta < 800MeV. If the elliptic flow originates 
from the collective flow of the fluid, the behavior at small momenta is linear 
in the transverse momentum (for particles with a small mass) 

'f2(P±) OCp_L . (3) 

Such a dependence is observed in hydrodynamic calculations (also including 
viscosity corrections) and in experiments on heavy-ion collisions at RHIC. 

The elliptic flow measured at the LHC is not increasing linearly with the 
momentum. The observed differential elliptic flow f2(p±) is larger than in hy- 
drodynamic calculations. The extrapolation of the measured elliptic flow to 
zero momentum is not going through the point ^2(0) = 0, as it should for the 
elliptic flow of a collective origin (Fig. H]). A linear flt to the data points with 
< IGeV gives an intercept of 0.06-0.01 ± 10~5, a small but statistically 
signiflcant and systematic deviation from zero. Moreover, the extrapolation 
of the results of hydrodynamic calculations would lead to a negative inter- 
cept. It is due to the fact that for massive particles we have f2(p±) oc 
While it is expected that hydrodynamic results deviate from the data at high 
momenta, it is very difficult to explain the excess of the observed elliptic 
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Fig. 3. Elliptic flow coefficient in Pb-Pb collisions as function of centrality, ALICE 
Collaboration data [Ij compared to viscous hydrodynamic results. 

flow at small momenta p± < SOOMeV. The data points lie above the lin- 
ear function from hydrodynamic calculations for all centralities in the range 
200MeV< p± < SOOMeV. The additional shift of the measurements above the 
calculated values is about 0.01-0.02. One possibility is that the linear behavior 
happens at very small momenta p < 200MeV, not observed experimentally. 
The data points observed by the ALICE Collaboration {pj_ > 200MeV) would 
correspond to the region of strong viscous correction, with deviations from 
the linear behavior. In the upper panel of Fig. H] is shown an example of a 
calculation starting with initial eccentricity increased by a factor 1.5 and with 
Tj/s = 0.24. Such a calculation gives the same integrated elliptic flow as ob- 
served experimentally. The calculated elliptical flow coefficient V2(p±) is closer 
to the data points, but still a systematic deviation is visible for small momenta. 

The elliptic flow in hydrodynamic models originates predominantly from the 
azimuthal asymmetry of the collective flow velocity. On general grounds, one 
expects a linear behavior of V2{p±) for light particles at small momenta. There 
is a different mechanism generating the azimuthal asymmetry, related to the 
asymmetry of surface of the emission. Due to strong jet quenching in the 
dense matter, the escaping partons originate from the surface of the flreball. 
In peripheral collisions the geometrical asymmetry of the source leads to pref- 
erential emission of escaping partons in-plane. The elliptic flow coefficient of 
the emitted partons would be approximately independent of the parton mo- 
mentum [16]. The final observed hadrons are distributed in a cone around the 
direction of the fragmenting jet parton. 

At the LHC energies, in each Pb-Pb collision many jets are formed. The jets 
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Fig. 4. Elliptic flow coeflicient of charged particles as function of transverse mo- 
mentum for centrality 30-40% (upper panel) and 10-20% (lower panel) measured 
by the ALICE Collaboration [1]. The dotted line represents viscous hydrodynamic 
results, the dashed line the approximate elliptic flow from jet fragmentation and 
the solid line the weighted mixture of the jet and hydrodynamic contributions (Eq. 
m. The dashed-dotted line in the upper panel represents the elliptic flow obtained 
increasing the Glauber model initial eccentricity by 50% and using rj/s = 0.24. 

interact with the dense matter in the fireball, leading to jet quenching and jet 
asymmetry [T7|[T8] . The interaction of the jet parton with the matter in the 
fireball, depends on the path it travels, which leads to an azimuthal asym- 
metry of the particle emission at large transverse momenta. The elliptic flow 
coefficient at large momenta at the RHIC energies is approximately indepen- 
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dent of the momentum and is of the order of 0.05-0.12, depending on the 
centrahty p^|20] . The value of the elhptic flow coefficients from jets at mo- 
menta p± < IGeV is not known experimentally, as this region is dominated 
by the statistical emission from the thermalized fluid. Once the direction of 
the jet parton is fixed, the fragmenting hadrons inherit, to a large extent, the 
original direction of the jet. A particle of momentum ~ 200-800MeV is not 
the leading particle in the fragmentation of the jet, it originates as one of the 
soft particles in the fragmentation of a parton with the momentum of a few 
GeVs. The leading particle is emitted in a narrow cone around the jet parton. 
The correlation between a high momentum trigger particle (leading particle) 
and associated particles show that also the soft particles are emitted in a cone 
around the jet direction. Although, one expects the jet cone opening is larger 
for particles with small momenta, we assume that the elliptic fiow coefficient 
from jets is approximately constant down to 200MeV for this first estimate. 

Following the parameterization of Ref . [2T] , we write the elliptic fiow coefficient 
at a given momentum as the sum of two components, a jet component V2^* 
and a hydrodynamic component ^2^'^^° 

V2{p±) = (1 - 9{P±))V2^'^'"{P±) + 9{P±)V2^\P±) , (4) 



g{p±) denotes the proportion of particles originating from jet fragmentation 
at a given momentum. It implicitly assumes that the jet particles, accounted 
for in the weight g{p±), are particles that do not thermalize and conserve their 
original (non- hydrodynamic) V2- There are only few constraints on the form of 
vl^* and g{p±)', g ^ 1 at large momenta and, from the deviation of the ALICE 
data away from the linear function V2{p±) oc p±, we have g{p±)v2'\p_i_) ^ 0.01- 
0.02 for p± < SOOMeV. For illustration, we fix the parameters [21j as g{p±) = 
(1 + tanh((p - p^)/Ap)) /2, p^ = 2.8GeV, Ap = 2.8GeV, VjetiPi) = 0.11 for 
centrality 30-40% and p^ = 3.2GeV, Ap = B.OGeV, Vjetip±) = 0.085 for cen- 
trality 10-20%. The final elliptic fiow coefficient from the two components is 
denoted by the solid lines in Fig. |H The deviation of the data from the hydro- 
dynamic calculation in the range 200MeV< p± < 800MeV can be explained as 
due to a 10-20% contribution of jet particles. For soft particles emerging from 
a jet parton, one expects some reduction of V2^'^{p±) for small momenta, due 
to a larger jet cone opening or to a possible coalescence with thermal partons. 
This would imply a larger value of g{p±) to get a similar V2'\p±)g{p±). 

We present results of relativistic viscous hydrodynamic calculations for Pb-Pb 
collisions at ^/s = 2.76TeV. The initial density for the evolution is fixed us- 
ing the experimental results for the charged particle multiplicity as function 
of the centrality. The same calculation was shown to describe correctly the 
observed interferometry radii in central collisions at the LHC [12]. The hy- 
drodynamic calculation using Glauber model initial profiles and a small shear 
viscosity coefficient reproduces the integrated elliptic fiow coefficient observed 
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by the ALICE Collaboration. Surprisingly, the differential elliptic flow as func- 
tion of the transverse momentum V2{p_l) shows systematic deviations from the 
hydrodynamics estimate. At small momenta the data points lie above the lin- 
ear hydro dynamic behavior V2{p±) oc p±. It is a very unexpected result. In 
the small transverse momentum region, besides the dominant statistical emis- 
sion from a collectively expanding fluid, a non- negligible contribution from 
jets appears. This observation, if conflrmed by further experimental studies, 
indicates that at the LHC energies non-thermalized particles from jet frag- 
mentation constitute a substantial part of the soft spectrum. This additional 
source of particles with soft momenta should be taken into account in quan- 
titative studies of hydrodynamic models of heavy-ion collisions, in calculating 
the spectra and the elliptic flow of particles. On the other hand, it could serve 
as a new frontier of studies of the jet parton formation, attenuation and frag- 
mentation in a heavy-ion environment. In particular, the fact that particles 
with soft momenta do not thermalize shows that the jet fragmentation occurs 
outside of the thermal fireball. 
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